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ABSTRACT 


The propagation of sound due to a line acoustic source in the 
moving stream across a semi-infinite vortex sheet which trails from a 
rigid plate is examined in a linear theory for the subsonic case. A 
solution for the transmitted sound field is obtained with the aid of 
multiple integral transforms and the Wiener-Hopf technique for both 
the steady-state (time-harmonic) and initial-value (impulsive source) 
situations. The contour of inverse transform and hence the decomposi- 
tion of the functions are determined through causality and radiation 
conditions. The solution obtained satisfies causality and the full 
Kutta conditions. The transmitted sound field is composed of two waves 
in both the steady-state and initial-value problems. One is the wave 
scattered from the edge of the plate which is associated with the bow 
wave and the instability wave. These bow waves and instability waves 
exist in the downstream sectors. The other is the wave transmitted 
through the vortex sheet which is also associated with the instability 
wave. This instability wave exists in a downstream sector, but if the 
line source is close to the rigid plate it is blocked by the plate and 
does not appear. The transmitted sound field can be divided into three 
regions. The first one is the region where the incident wave is shaded 
by the plate and only the wave scattered from the plate edge exists. 

In the second region there exist the waves transmitted through the 



vortex sheet in addition to the scattered waves from the edge of the 
plate. The third region is the transition region between the two 
regions just described. In this region the waves transmitted through 
the vortex sheet can be heard partially depending on its wave number. 
The asymptotic nature of this region is the same as that of the second 


region. 
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CHAPTER 1 


INTRODUCTION 


The influence of a moving stream on the propagation of sound is a 
problem of some practical importance and has attracted a number of 
investigations in recent years. For the sake of analytical simplicity, 
most theoretical investigations related to this problem have been 
confined to the model study of two infinite two-dimensional, inviscid 
fluid half-spaces in relative motion separated by a plane vortex sheet. 
Sound waves emitted from impulsive (harmonic) line or point sources 
will impinge on the vortex sheet and be reflected and transmitted to 
the two half-spaces as reflected and transmitted sound fields. 

For linearized theory, multiple-Fourier transform techniques are 
most conveniently used to obtain the analytical solutions for these 
sound fields. While it is almost trivial to find the transformed solu- 
tions, the complexities of inverse transforms present much difficulty 
in finding a solution in a form nx>st suitable for practical applications. 

Jones and Morgan and their associates [ 1—5] in a series of papers 
investigated the problem of line and point sources. As they carried 
out inverse transform with respect to wave number first along the real 
axis, their solutions did not satisfy causality. To satisfy causality, 
they had to add the homogeneous solution which corresponds to instability 
waves. Chao [6] has indicated that if the inverse transform is carried 
out with respect to frequency first, the proper deformation of contour 
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which brings the causal solution is obvious and the instability waves 
are automatically included in the solution obtained. 

In the real world a vortex sheet between two fluids in relative 
motion can only exist downstream of a plate which separates them. The 
problem becomes a mixed boundary value problem and the Wiener-Hopf 
technique [7-8] has to be used to find the transformed solution. 

Orszag and Crow [9] have investigated this problem in which there are 
no external sound sources and the flow is incompressible . Cnghton [10] 
has extended their analysis to compressible subsonic flow. Morgan [ll] 
and Crighton and Leppington [12] have investigated the reflected sound 
field produced by a line source near a semi-infinite vortex sheet which 
trails from a rigid plate where a line source is situated in the still 
fluid. 

All of these problems require the decomposition of a known function 
into plus and minus functions which are analytic in the lower and upper 
wave number plane, respectively. However, due to the fact that there 
are poles whose location is frequency-dependent in such a way that it 
will move from lower half-plane to upper half wave number plane according 
to the value of frequency, the determination of a common region of 
analyticity and hence the contour of the inverse transform becomes a 
rather complicated matter and consequently requires additional assumptions 
to ascertain the desired contour of the inverse transform. Morgan has 
enforced the condition that the slope as well as the displacement of the 
vortex sheet are assumed to be zero at the plate edge (full Kutta condition) . 
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Crighton and Leppington have decomposed the function by assuming a 
large imaginary part for the frequency. They both have shown that 
their assumptions require a deformation of the contour of the inverse 
transform and satisfy the criterion for causality established by Jones 
and Morgan [13]. This criterion requires the analyticity of the solution 
in the upper half frequency plane for causality and has been used suc- 
cessfully by Morgan [14] and Munt [15] in their analysis of the cylin- 
drical vortex sheet problem. 

In Chapter 4 it is shown that the deformation of the contour, and 
hence the common region of analycity can be determined through causality, 
exponential decay of the functions , and radiation condition. In all 
three of these approaches, the decomposition of the function is unique. 

In Chapter 5 we determine the value of the entire function, 
assuming the edge condition that the component of the disturbed velocity 
in the stream direction is finite at the edge of the plate. This 
assumption requires that the full Kutta condition be satisfied as a 
result. Howe [16] has found that for a very low Mach number flow case, 
the satisfaction of the Kutta condition requires that external forces 
be present m the vicinity of the plate edge. Similar situations 
certainly would occur for the present problem. The existence of this 
external force probably can be traced to the viscous effect which should 
not be neglected due to the large velocity gradient in the vicinity of 
the plate edge. A detailed study of viscous mixing at a trailing edge 
by Daniels [17] suggests that the full Kutta condition is the logical 
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condition to apply at the plate edge. A further justification of the 
Kutta condition may be found in the experimental work by Bechert and 
Pfizenmaier [18]. Their work suggests that the full Kutta condition is 
the correct condition to apply if the shear layer is thin, which is the 
case in our problem. 

In Chapter 6 the solution is compared with that of the Sommerfeld 
half-plane diffraction problem (without the flow) . It is found that 
the existence of the flow smooths the edge behavior of the solution. 

In Chapter 7 the steady-state solution is calculated. The trans- 
mitted sound field can be divided into three regions depending on how 
the waves transmitted through the vortex sheet are shaded by the rigid 
plate. In Region 1, the incident wave is shaded completely by the 
plate and only the waves scattered from the plate edge exist. In 
Region 2 the waves transmitted through the vortex sheet exist m addition 
to the edge scattered waves. Region 3 is the transition region between 
the two where the part of the waves transmitted through the vortex sheet 
exist depending on its wave number. The existence of Region 3 is the 
direct conclusion of the fact that the line source and observer are 
located in the different mediums. It is expected that a similar situation 
will occur if the speed of sound is different between the two mediums 
even though the velocities of the flows in the two mediums are the same. 

In Chapter 8 the asymptotic evaluation of the transmitted sound 
field is carried out as far as possible. In the case where the line 
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source is situated close to the plate edge in comparison to the observer, 
the asymptotic nature of Region 3 is similar to that of Region 2, and in 
the far field the boundary between Regions 2 and 3 vanishes. 

In Chapter 9, the solution for the initial value problem is 
calculated. The wave fronts of the diffracted waves are obtained 
explicitly. 
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CHAPTER 2 


FORMULATION OF THE PROBLEM 


Consider the problem m which a moving fluid m the half-space y<0 

is separated by a rigid plate and a vortex sheet from a still medium m 

y>0. The rigid plate occupies the half -plane y=0, x<0, the vortex 

sheet occupies y= 0, x>0 in a Cartesian coordinate system and it is 

assumed that the flow is in the direction of increasing x. A line source 

of unit strength is situated in parallel to the x-axis through (x^y^) 

in the moving medium, so that y < 0 (Figure 2-1) • For the sake of 

o 

simplicity, the sound speed and density are assumed to be the same m 
both fluids- Because of the symmetry of the problem all quantities are 
independent of z. 


y 
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Figure 2-1 Sketch of the Model Geometry 
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Let and <t> 2 denote the velocity potentials of the small disturbed 
motions m the two fluid spaces. Due to the excitation of the line 
source , the governing differential equations for <f>^ and <J> 2 ar e: 

2 

(it + v ik) *2 " <t> 2 = <S(x_x o ) 6 (y-y 0 ) 6 (y < o) 

a 2 <{) 

— - avl = 0 (y> 0) 

9t 1 

where a is the speed of the sound and V is the velocity of the moving flow. 



Define the Fourier transform of . as: 

l 


f f'j**-*^ 


-i(kx-wt) 


dxdt 


( 2 - 2 ) 


Then Eqs. (2-1) reduce to the form: 


-lkx 


a 4> 2 2 = -5(y-y o )e 


+ = 


. 2 ’ 1 2 T 2 .2 2 

3y 4TI a 


a 2 ^ 


3y 


2 l+ ^l = ° 


1 

1 


(2-3) 


where 
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and M = — is the Mach number of the lower medium. Let y = h(x,t) 
a 

indicate the equation of the vortex sheet. The boundary conditions on 
the vortex sheet are the continuity of pressure, giving 





9^ 

3t~ ' 


(2-4) 


and the kinematic condition of the continuity of particle displacement 
which implies 


j 3 

/ 3<) V 

3 ^2 

3 ) 

| 3t + 

( V + 3x , 

I 3x + 3y 

3y j 


H(x,y, t) = 0 , 


(2-5) 


where H(x,y,t) = y-h(x,t). Linearizing the equation, we have 
34>. 


/ 3 3 3 \ 

(^ +V ^ + 9T37) (Y " h) 


= 0 


and 


3h _ 3h ^^2 3y 

3t 3x 3y 3y 


= 0 


We get: 
3y 


9h + 

3t 3x 


(2-5a) 


Similarly, for medium (I) 



3y 3t * 


(2— 5b) 
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The boundary condition for x < 0 is the vanishing of the normal velocity 


on the splitter plate which implies 


3 ^ 3<(>2 

3y ~ 3y 


= 0, h = 0 . 


Define the half -range Fourier transforms of ^ and h as: 


‘f’i+OWYfW) = -ij j* J ^(XrY/tJe 1 


-i (kx-cot) 


dxdt 


-oo 0 


> + (k,w) = J * J “ h(x,t)e i 


-i (kx-wt) 


dxdt 


-oo 0 


( 2 - 6 ) 


(2-7a) 


(2-7b) 





OO 0 

J’ J ' (^(XfYjtJe 1 ^ kx "° Jt ^dxdt 

— oo — OO 


(2-7c) 


a_(k,u>) = -^2 j r J ” h(x,t)e 1 


-l (kx-wt) 


dxdt 


(2-7d) 


Equations (2-4) , (2-5a) , (2-5b) , and (2-6) reduce to the form 


- 1W V + 2? 


OO ^ 

9<J> 


f 3 *2 

J (x ' y ' 


to) e 1 ^ cx dx = -ioxJ>. 


1 + 


Vi 


2+ 27T03 


^ + 


-r -ikx 

^2 e 


+ ik 


OO 

/ 


-r -ikx, 

4>2 e dx 


"*1 + ' 


^2+ + L 1 



(2-8) 
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where 


, Vk 
a = 1 -- 


L i“’isrV x ' y '“ ) 


x f y=0 


dtp 


2 + 

dy 




3h (x,y,a)) -ikx, 
e d * 



r 00 ^ l 

- - * 57 

— -ikx 
h e 

+ lk / he ikx dx 

J 

C 

) 0 -I 


itoah + . 


(2-9) 


Here, we use the condition h| =h| = 0, because at the edge of the 

x=0 x=0 

plate , the vortex sheet must attach to the plate. Similarly, 


3<J> 


1 + 


3y 


= - io)h 


( 2 - 10 ) 


and 


3<J> 


1- 


3y 


3<J) 


2 - 


3y 


= 0 h = 0 


( 2 - 11 ) 


Solutions to the set of the ordinary differential Eq. (2-3) can be 
easily found as : 


iY x y 


<p^ = A 1 (k / to)e + B 1 (k,a))e 


-iY L y 


-ikx 

lY 2 y e ° 

4> 2 = A 2 (k,to)e +B 2 (k,aj)e — ^ — sin 

4ir a Y- 


-iy 2 y 


in [ Y 2 (y_y o ) ] 


(2-12) 


H(y-y ) 
o 
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where H( ) is the Heaviside step function. The sign conventions of 
and Y 2 are chosen such that the signs of their imaginary part are always 
positive. Consequently, and terms of the solution are finite at 

infinity for all values of k and to. and must be set to zero. 

Using the relation 


3y 

y=0 


iY A I 

y=0 


3y 

y=0 


= -iY <$> 
t 2 y 2 


y=0 


-i(kx +Y o y ) 

O 2 0 
e 


. 2 2 

4ir a 


d> = <b. + <f>. 

i+ i- 


and the boundary conditions (2-9) , (2-10) , and (2-11) , we obtain 



y=0 



y=0 



y=0 




y=0 


-i(kx +Y y ) 

O 2 O 
2~2 

4tt a y 2 


(2-13) 


Substituting these into Eq. (2-8) gives the Wiener-Hopf equation as: 


(Dh + = K ^ 2 _ 


- <j> 


1- 


y=0 y=0 




iy.a e 


-i(kx o+ Y 2 y o ) 


2 2 2 
4ir a (Yj+Y^ ) 


(2-14) 


where 


K = 


y 2 

W 2 


ii 



If we can decompose K into plus and minus functions as: 


K 



we can rewrite the Wiener-Hopf equation as: 


a)h + K + = K_ 


-♦i- 

y=0 y=0 


- L i'- 


-i (kx +y y ) 

. o z o 

iYjOl e K + 

2 2 2 
4ir a (Y^Yja ) 


(2-15) 


The last term of Eq. (2-15) can be decomposed into the sum of plus and 
minus functions as: 


F 


-1 (kx +Y,y ) 

O Z O 

lYja e K + 

2 2 2 
4tt a (Y 2 +Y 1 a ) 


= F + F 
+ 


where F + , F are given as: 


F + (k) 



F(X) 

X-k 


dX 


F (k) 



F (X) 
X-k 


dX . 


(2-16) 


I (2-17) 


The exponential decay of F(k) on the contour, since Inry >0, y <0 

z o 

ensures that these integrals exist and that F + (k) =0(|k| , F_(k) =0(|k| ^ ) 

as | k | ■* 00 - The meaning of the plus and minus functions and the decompo- 
sition will be discussed in the next chapter. 
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Substituting Eq. (2-16) into Eq. (2-15) completes the decomposi- 


tion, to give: 


(Dh + K + + F + = K_ ^ 2 _ 


- <t> 


i- 


y=0 y=0 


- L, " P, 


= C(k) . 


(2-18) 


The function C(k) defined by Eq. (2-18) is an entire function of k and 
must be a regular function of k in the whole k-plane. If C(k) is 
determined, 0)h + is given as: 


wh + = 5T ( c(k) _F +) • 


(2-19) 


From Eq. (2-13) 


-C(k) + 
Al = Y 1 K + 


-ikx 


( 2 - 20 ) 




2 2 . COSY 2 y o 


4tt a Y 


The solutions in the transformed region are determined. 
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CHAPTER 3 


DECOMPOSITION OF THE FUNCTIONS 


First, let us construct the branch cuts of and in the k-plane. 
To do this, the sign of ImO) has to be investigated. One of the ways to 
determine the sign of Imto is to consider the wave equation with 
damping [8]. 



3 2 <j> 

at 2 


£ 34 > _ 

2 at" 


o , 


e >o . 


(3-1) 


Carrying out the Fourier transform, defined as: 
00 

* = 2rr J * e dt ' 

— oo 

Eq. (3-1) reduces to: 


V (}) + " <P = 0 , 


(3-2) 


(3-3) 


where (0 is assumed to be real. Alternatively, assuming to to be complex 
as (0 = ( 0 ^ + ito^ , where to^ and to 2 are real and taking the wave equation 
which does not have damping, 



a 



(3-4) 


14 



Carrying out the Fourier transform, we obtain 


2 2 „ 

0 ) - 0 ) + 203 0)1 

+ — 2 2 1 2 <J> = 0 

a 


Comparing Eq. (3-3) and (3-5) , it is concluded that assuming to is real 
and considering the wave equation which includes damping is equivalent 
to assuming that 0) is complex and considering the wave equation without 
damping, if lma)>0. Therefore we should assume lmo>0. 


The sign of Imo) can also be determined through the causality 
condition which requires that no disturbance should be present before 
time t = 0. In order to satisfy this, the contour of inverse transform 
has to lie in the upper-half domain of the CD-plane, which indicates that 
lma)>0. The causality condition will be discussed m Chapter 4. 


Consideration of the radiation condition gives another way to 
determine the sign of ImO). Consider the branch cuts of y in the k -plane. 
For simplicity, assume Re(i)>0. Inty^ must be positive according to the 
commitment made in Chapter 2, Two kinds of branch cuts are possible as 
shown in Figure 3-1. For the case ImO)>0, Rey^ is positive on the real 
axis and considering the form of the inverse transform: 


- II 


iy^+ikx-iwt 

U))e da)dk 


(3-6) 


this indicates the outgoing wave in y-direction. On the other hand, for 
the case lmu)<0, Rey^ is negative on the real axis, which indicates the 
incoming wave in y-direction. Therefore Inti) should be positive. 
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indicates the argument of ^^-+k 
□ indicates the argument of k 

Figure 3-1 Branch Cuts of in k-Plane 


In all these methods, the common idea is to set a direction in time 
space, from the past to the future. 

The poles of K are discussed in detail in Appendix A. The results 
are as follows: When 2 >M, there are two complex poles at: 


k = u(M)k 


and 


k = u(M)k =u* (M)k 
o o 


where u(M) indicates the complex conjugate of u(M) and 

+ Vm 2 +4-4 Vm^H 


u(M) = 


M 


+1 


:V » 2 


2 VM +1 - 2 


(3-7) 

(3-8) 


(3-9) 


k = to/a . 
o 
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The poles k = u(M)k , k=u*(M)k are associated with instability of the 
o o 

vortex sheet when M < 2/2 and play a prominent role in the solution of 
the problem. They are therefore displayed explicitly by defining y(k) by 


K(k) 


Y 1 (k) Y 2 (k) 

y(k) (k-u(M)k ) (k-u*(M)k ) 

o o 


(3-10) 


We seek the radiating acoustic solutions that decay like 

-Imk | x | 
o 
e 

or 

-Imk | x | / (1 +.M) 
e 


as x + + oo. Thus, all plus functions, denoted by a subscript + , will be 
analytic in the region 

Imk 

Imk < T+m ' (3_11) 


while all minus functions, denoted by a subscript -, are analytic for k 
in the region 

Imk > - Imk . (3-12) 

o 

There is evidently a strip 
Imk 

\ > Imk > - Imk : Region A (3-13) 

1+M o 

in which full-range transforms of (f>^, (j> 2 , and h are analytic functions 
of k. 
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We consider the case M<1. In this problem the branch cuts of 

and should take the form as shown m Figure 3-2 to satisfy the 

commitment that Imy^ and Imy 2 must be positive on the real axis. The 

branch cuts of /k -Mk+k and Sk +k lie m Region B and the branch cuts 
o o 

of /k -Mk-k and /k -k lie in Region C even though U) changes it value as 
o o 

long as lmu)>0, where Region B and Region C are defined as: 


Region B: link < - Imk^ 

link 

Region C: Imk > 7 * 

1+M 



Region C 


Region A 


Region B 


Figure 3-2 Branch Cuts and the Common Region of Analyticity 
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Figure 3-3 


The Position of the Poles m the k-Plane 
According to the Value in the co-Plane 
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The functions y(k) , Y , Y are split into two parts such that: 

l ^ 

U(M - u + oo/u_(M , Y 2 =Y 2 ./r 2+ . 


From the branch cuts we made previously it is apparent that 


v (k) 


<k -k> 
o 


- 1/2 


Yj_(k> 


Y 2t (W 


(k +k) 
o 


1/2 


(k -Mk -k) 
o 


- 1/2 


(3-14) 


Y 2 Jk) 


(k -Mk +k) 
o 


+ k> 1/2 


The calculation of y + (k) and y (k) is presented in Appendix B. The 
difficulty arises from the poles k = u(M)k Q and k = u*(M)k Q . As can be 
seen in Figure 3-3 , the position of the poles in the k-plane changes 
according to the value of 0). As a result, we cannot determine whether 
the terms £k-u(M)k Q J, £k-u*(M)k o J work as a plus function 
function. 


or a minus 
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CHAPTER 4 


CAUSALITY AND RADIATION CONDITION 

The causality condition requires that no disturbance be present 

before time t = 0. Let the transformed solution have two poles at 

k=u(M)k and k = u*(M)k , which is the case in our problem. First, 
o o 

consider the inverse transform with respect to to, assuming that k is 
real. It is apparent that to satisfy causality, the contour must lie 
above all the singularities in the to-plane. In other words, the contour 
is deformed so that it always includes poles as shown in Figure 4-1. 

Then the branch cuts are constructed so that the radiation condition 
is satisifed along its contour; namely, 

Iirry^ > 0 along the real axis and the position of the 

pole enclosed by the deformed contour, (4-la) 

Recd-Rey^ >0 along the real axis and the position of the 

pole enclosed by the deformed contour. (4-lb) 



Figure 4-1 Contour of Integration and Branch Cuts in a)-Plane 
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The condition (4-lb) indicates that the solution leads to outgoing waves 
in y >0. It is not easy to express the radiation condition for y<0 
because of the moving fluid. However , one can determine the branch 
cuts of in such a way that when M^O they coincide with those of y . 
These branch cuts are also shown in Figure 4-1. To obtain the final 
result, integration is carried out with respect to k along the real axis. 


Next, consider the inverse transform with respect to k, first 
assuming that 0) is real or has a very small positive imaginary part. 


f (x,03) 



f (k,0})e^ X dk 


(4-2) 


If we integrate it along the real axis, f(x,oj) has singularities in the 

03-plane when the poles k=u(M)k , k = u*(M)k hit the contour as 0) changes 

o o 

its value as can be seen in Figure 4-2. 



Figure 4-2 Movement of the Poles as 03 Changes from -<» to » 
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There is no proper contour m the 03-plane which satisfies causality because 
of these singularities. To overcome this, the contour C is deformed in 

JC 

such a way that after the pole hits the real axis in the k-plane, the 
contour always includes this pole as shown in Figure 4-3 and as a 
result the singularities disappear. The branch cuts in the k-plane 
have already been made in Chapter 3 and along the real axis the 
conditions (4-1) are satisifed. We have no choice in making the branch 
cuts to satisfy conditions (4-1) at the position of the pole enclosed 
by the deformed contour except for the geometrical relation between 
the branch cuts and the poles. However, we have a choice in the defor- 
mation of the contour. The contour may be deformed to include the pole 
in the upper-half plane as shown in Figure 4-3 (a), or the contour may 
be deformed to include the pole in the lower-half plane as shown m 
Figure 4-3 (b). 



(a) 


(b) 


Figure 4-3 Contour of Integration in k-Plane when Reoi > 0 
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Consider the case Reu)>0. The relation between the two branch 
cuts of y and two poles is as Case I or Case II shown in Figure 4-4 (a). 
According to the branch cuts defined in Chapter 3 at the pole k = u(M)k^ 
we have : 


Case I 


arg(k 

o 

-w i/2 

+ 

II 

with 

i> 

2 

0 > 0 
a 

(4-3a) 

arg(k 

o 

+ M 1/2 

II 

H “ 

with 

I> 

2 

0. > 0 
b 

(4-3b) 


and 


7T+0 +0, 
a b 

arg Y 1 = 2 


7T > arg y 1 > - , 


(4- 3c) 
(4-3d) 


Case II 


arg(k -k) 
o 


1/2 


0 

tt a 
2 + 2 


ar,( k + W 1/2 .% 

O 2 


and 


arg Y 1 = 


— TT+0 +0^ 
a b 


0 > arg Y 1 > - j . 


(4-4a) 

(4-4b) 

(4-4c) 

(4-4d) 


Therefore at the pole k=u(M)k , condition (4-la) is satisifed but (4-lb) 

o 

is not satisfied in Case 1, and condition (4-lb) is satisifed but (4-la) 
is not satisfied in Case II. 
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(a) Rem > 0 




(b) Red) < 0 


Figure 4-4 Relation Between Branch Cuts and Poles 
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At the pole k = u(M)*k , both Case I and Case II give: 

o 


arg (k 

o 


arg(k 

o 



+ k) 


1/2 



and 


7 T — 0 -0 
a 


b 


arg Y 1 = 



arg 


2 

> 0 . 


(4-5a) 

(4-5b) 

(4-5c) 

(4-5d) 


Therefore, at the pole k = u*(M)k , 

o 

and (4-lb) • 


satisfies both conditions 


(4-la) 


Next, consider the case ReO) < 0. The relation between the branch 
cuts and two poles is as shown m Figure 4-4 (b) . Similarly, according 
to the branch cuts defined in Chapter 3 at the pole k = u*(M)k Q we have: 


Case I 

arg(k -k) 1/2 = 
o 

arg(k +k) 1/2 = 
o 

and 


arg Y x = — 
ir . 

2 > ar 9 Y l 


-0 
c 

2 

> 0 



with 1 > e„ > o 

2 c 

with ~n > > 1 


(4-6a) 

(4-6b) 

(4-6c) 

(4-6d) 


26 



Cage II 


arg(k -k ) 1/2 = tt - -f 

o Z 

1/2 ®d 
arg(k + k) i/Z = — 
o z 

and 

27T-e +6, 

c a 

arg 5 

it > arg y 1 > it . 


(4-7a) 

(4-7b) 

(4-7c) 
(4— 7d) 


At the pole k=u*(M)k , condition (4-la) is satisfied but (4-lb) is not 

o 

satisfied in Case I, and condition (4-lb) is satisifed but (4-la) is not 
satisifed m Case II. 


At the pole k=u(M)k , both Case I and Case II give: 

o 


arg (k -k ) 1/2 = ~ 
o 2 


arg (k +k) 1/2 = 7 r - 

o z 


and 


27T+0 -0 
c d 


arg y ± - - 2 


it > arg Y 1 > j • 


(4-8a) 

(4-8b) 

(4-8c) 

(4-8d) 


Therefore, at the pole k = u(M)k , y satisfies both conditions (4-la) 


and (4-lb) . 
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It is concluded that if we deform the contour so that all the 
poles are above the contour, radiation conditions (4-1) are satisfied at 
the position of the poles k=u(M)k Q , k = u*(M)k o , as well as along the 
real axis, whereas if we deform the contour so that all the poles are 
below the contour, radiation conditions are not satisfied. Therefore, 
the contour must be deformed to include all the poles above it. 

It is obvious that the contour obtained here is also valid for the 
infinite vortex sheet problem and it is easily shown that the contribu- 
tion from these poles coincides with the homogeneous solution obtained 
by Jones et al. [1-5]. 

The decomposition of the functions in the Wiener-Hopf technique is 
slightly changed. As the contour of inverse transform must lie m the 
common region of analyticity, the common region of analyticity has to be 
deformed m a similar way. Plus functions are defined as the functions 
which are analytic below the contour and minus functions are defined as 
the functions which are ana lytic above the contour. Then the decompo- 
sition of K(k) is determined uniquely as follows: 

K + = Y 1+ Y 2+ y + (k) |k-u(M)k Q | jk-u*(M)k Q j- (4-9a) 

K_ = y 1 _ y 2 _ y_(k) . (4-9b) 


28 



CHAPTER 5 


THE CONDITIONS AT THE PLATE EDGE AND THE VALUE OF THE ENTIRE FUNCTION 


The Wiener-Hopf equation is given in Chapter 2 as: 


Wh + K + + F + = K_ (a<t> 2 _ 


-* 1 - 
y— 0 y=0 


- L i)- r - - 


C(k) 


(5-1) 


where 




L 1 27T0) < f ) 2 (X/CA)) 


x=0,y=0 


C(k) is the entire function of k and plus functions and minus functions 
are analytic below and above the contour determined in Chapter 4, 
respectively. Near the edge of the plate at the origin .we assume: 


$^(x,U3) C^(U3) as x -0 on y = +0 (5-2a) 

$ 2 (x,0)) C 2 (w) as x -0 on y = -0 (5-2b) 

— — 

(J> 2 (x,0)) - $ 2 ( 0 , 03 ) = 0 (x ) as x + +0 on y = -0 (5-2c) 

— ^2 

h(x,03) = 0 (x ) as x +0 on y = +0 (5-2d) 


where £^>_1, ^>0. In these conditions, the expression x^-0, y = +0, 
for example, means that x tends to zero through negative value of x on 
the upper side of the plate. The C^(U3) are functions of 0) and need not 
be known explicitly. 
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The third condition indicates that the x-component of the disturbed 
velocity is finite below the vortex sheet at the edge. The fourth 
condition is the requirement that the vortex sheet should attach to the 
edge of the plate. 


Carrying out the half-range Fourier transform defined by Eq. (2-7) , 
the asymptotic behavior of the above functions can be calculated with 
the aid of the Abelian theorem as: 


4> 1 _ (k,(d) = o (|kl -1 ) 

7 2 _0c,u>) = 0 (|k| -1 ) 


as k -*•«* in Imk > -link 


as k • >co in Imk > -Imk 


(5-3a) 
(5- 3b) 


= i<M0,W) 

*2- (k ' W) ’ 27Tk = ° l|k 


h + (k,w) = 0 | |k | 2 j 


(|ki L j as k 


in Imk > -Itfik 


as k m Imk < 


Imk 
o 

1+M 


(5-3c) 


(5-3d) 


F + (k) , F (k) are given as follows: 


? + (k > ’Si / T5T dX 


(5-4a) 


F (k) 


--L- f 

2tti J 


FO) 

X-k 


dX 


(5-4b) 


where 


F (X) = 


^ (X)cte 


-i(Xx+y(X)y) -i(Xx +Y (X)y ) 

' O 2 O / _ _ 2 /*v\ ' O 2 O / 

K + iy_(X)Y 1 _(X)Y 2+ (X)ae 


4ir 2 a 2 (y 2 ^) + Y 1 (^)« 2 ) 


.2 2 

4ir a 


(5-4c) 
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Since F(X) does not have a pole at X = u(M)k nor X = u*(M)k , the contour 

o o 

C + and C_ need not be deformed and lie in Region A defined m Chapter 3. 

C + passes above X=k while C passes below X = k. The exponential decay 

of F ( X) on the contour, since Imy. > 0, y <0 ensures that these 

z o 

integrals exist and that F + (k) = 0 (|k| F_(k) = 0 ( |k | as k in 

respective regions. The asymptotic behavior of y + (k) , u (k) is calculated 
in Appendix B as: 

y + (k) = 0 (|k | 1/2 ) 
y_(k) = o(|kf 1/2 ) 


Imk 

as k -*-°° in Imk < , ° (5-5a) 

1-rM 

as k “►«> m Imk > - Imk . (5-5b) 


Therefore, the asymptotic behavior of the known functions in the Wiener- 
Hopf equation can be estimated as follows: 


K + (k) = 0 ( 

>1 

1 3/2) 

as 

k 

K_(k) = o( 

l k 

,l / 2) 

as 

k ■><» 

o 

II 

a 

+ 

Pm 

l k 

r 1 ) 

as 

k -*<» 

F_(k) = 0 ( 

> 

i- 1 ) 

as 

k ->co 


Imk 


in 

Imk 

< 

o 

1+M 

(5-6a) 

in 

Imk 

> 

- Imk 

(5-6b) 




o 





Imk 


in 

Imk 

< 

o 

l+M 

(5-6c) 

in 

Imk 

> 

- Imk 

(5-6d) 


o 
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First we calculate the value of minus functions: 


L i - - 1ST *2 (0 -“> 


- - T V k '“ | l +0 

y =0 


(n" 1 ) 


(5-7) 


c_ (a? 2 _ -4> 1 _ -L 1 ) - F_ =K_ U 2 _ +o(|k| X jJ -F_ 


y =0 y =0 


y =0 y =0 


/ 1/2 "^\ 

= 0 ^|k| J , as k 


->-00 in imk > - Imk 


(5-8) 


As l >_ 1, the value of the entire function C(k) must be identically 
zero because of Liouville's theorem. Next, calculate plus functions 


l K + F = 0 ( Ik 

+ + + y 


\/ 2 ~%2 \ Imk Q 

) as k -+ OT in Imk < - 77 — . 
/ 1+M 


(5-9) 


In order to make C (k) equal zero, > l/ 2 * If we set C(k) =0, then 


— r 

1 = — 7 - = 0 ( | k | 5//2 ) as k -+ 00 m Imk < 

t K . ' ' 


1+M ' 


(5-10) 


which implies 


h (x,oj) « o(x 3//2 j , 


as x -*■ +0 on y = 0 , 


(5-11) 


and as a result the slope of the vortex sheet becomes zero at the edge 
of the plate as well as the displacement of the vortex sheet (full Kutta 


condition) . 
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CHAPTER 6 


BEHAVIOR OF THE SOLUTION WHEN M TENDS TO ZERO 


The behavior of the solution near the edge of the plate is different 
from that of the Somxnerfeld half-plane diffraction problem (without the 
flow) . Let us examine the y-component of the velocity on the vortex 
sheet near the edge of the plate. 


K + - Y 1+ Y 2+ P + 00 - u(Mlk Q J- |k-u*(H}k o | 

’ V2 + ( Y 2 + V 2 ) U - lkl ' 

since 


1 

(w 2 

)y_oo 

I* 

-u(M)k Q j 

|k-u* (M)k Q j 


( 6 - 1 ) 


(6-2) 


The asymptotic behavior of y_ (k) is given in Appendix B as : 

y_(k) = o(m - 1 |k 1 -1/2 ) ; (6-3) 

Hence 


K + (k) = 0 | M -1 |k I 1/2 (l+M|k|+M 2 |k| 2 )| . (6-4) 

The asymptotic behavior of F + (k) is given in Chapter 5 as: 

F + (k) = 0 (m"* 1 | k | _1 ) . (6-5) 
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Therefore in our problem 



(6-6a) 


(6-6b) 


(6-7a) 


(6-7b) 


Here the dependence of the asymptotic behavior on M is displayed 
explicitly. If M is finite 



(6-8a) 



(6-8b) 
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Then the Abelian theorem implies 



as x++0 on y = +0 (6-9a) 


as x + +0 on y = -0 . (6-9b) 


In the Sommerfeld problem 



as x +0 on y = 0 . 


(6-10) 


( 6 - 11 ) 


In these expressions the asymptotic behavior is calculated when k*><» in 
the lower half-plane of the k-plane defined in each problem respectively. 
If M+0 and Mk+0, the asymptotic behavior of the two problems coincides. 
In other words, the existence of the flow smooths the edge behavior of 
the solution. In fact, for the small Mach number 


u(M) « 

M 


(6-12a) 

u* (M)ss^ 
M 


(6-12b) 

Y^(k) «Y 2 (k) «k 2 -k 2 


(6-13a) 

Yl + 0O « y 2+ (k) « (k Q -k] 

| — 1/2 

(6-13b) 

Y 1 _(k) « Y 2 _(k) ~ (k o+ k) 

i 1/2 

(6-13c) 
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y (k) 


(V k ) 


- 1/2 


M 


(V k ) 


1/2 


M 


U. (k) w 2 

k 


(6-13d) 


(6-13e) 


Hence 


M 


K + (k) = 


(k-u(M)k o J |k-u*(M)k o | 


? (k) = - f 

T / 8tt a 


(k -k) 1/2 
o \ o / 

- i ( X V Y 2 (X)y o) dX 


(6-14a) 


3 a 2 (k Q -X) 1/2 (X-k) 


C+ 


4> 1 (x,y,a)) = 


(6-14b) 


M 


-l (Xx q +y 2 (XJy^kx-^ (k)yj 


1 1 O” 

J O 3 2, 

</ 8tt a | 

* 

o 

1 

>-* 

| 1/2 <A-MH 2 | 

|k-u(M)k Q j 

M 

[k-u*(M)k ) 
^ of 

<V k ) 

l 1/2 


(6-14c) 


As M 0 and Mk -> 0 , 
a +1 

and 

M 2 (k-uk ) (k-u*k ) - 

^ 2-« 1 +a 2 ■+ 2 


(6-15a) 


(6-15b) 
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We have 


c k c + 


-i (Xx q +Y 1 (X) y Q -kx-Y ;L (k) yj 


dXdk . (6-16) 


* 

The solution coincides with that of the Sommerfeld problem discussed in 
Appendix C. 
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CHAPTER 7 


TRANSMITTED SOUND FIELD 


The transmitted sound field is formally given by the equation 


*1 



C 

0) 


F + (k) 

VT 


ikx-icut+iy y 
e 1 dkda) 


( 7 - 1 ) 


where and are the contours discussed m Chapter 4 and 


V k > -sr/iSr di <7 - 2a > 

c + 

ly (X)y, (A.) y (X)a(X)e 

F(X) = — — - — (7-2b) 

4tt a 


where C + passes above X=k. If we change the contour of the integration 
so that C + passes below X=k, the contribution from the pole X=k must 
be added to the integration, which gives 


F + (k) = F (k) 



F (X) 
X-k 


dX 


= F (k) - F (k) . 


(7-3) 
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The transmitted sound field 0^ due to F(k) is given as 


•.-// 


lY l- Y 2+ P - (k)ae 


ikx-iajt+iy v-ikx -iy_y 
1 o 2 o 


C C, 
a) k 


2 2 

4tt a Y 1 K + 


dkdoj 


// 

C co C k 


iae 


lkx- ia)t +i y^y- ikx^ - i y 2 Y c 


2 2 2 
4tt a (Y 2 +Y 1 a ) 


dkda) , 


(7-4) 


which is exactly the same as the solution for a doubly infinite vortex 
sheet problem without a rigid plate as can be shown in Appendix D. 
Therefore, F (k) is considered to indicate the effect of a semi-infinite 
plate on the doubly infinite vortex sheet sound field . 


In this chapter we carry out the integration (7-1) with respect 


to k. The integration with respect to U) is carried out in Chapter 9. 


The advantages of integrating with respect to k first are: 


(l) After the integration with respect to k, the time harmonic 

solution is obtained. The time harmonic factor is assumed 

. , -icot 

to be e 


(n) The geometrical relation between the source and the 

observer is more clearly investigated as will be shown 
later. 


Introducing the new parameters 

C = kA 0 

V = A/k , 
o 


(7-5) 
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the time harmonic solution of the transmitted field can be written as : 


4 ) 1 (x,y,o)) = 


v ^(V)w 1 _(v) w 2+ (V) (l-MV)e 


/• /• v _(v)w 1 

J J 87T 


-ik Vx +w (v)y -£x-w_ (?)y 
002 o 1 


c C C v 


8lt a (V-C) (?-u(M)) (c-u*(M))w 1 _(C)w 2+ (?)v + (0 


where w.^ (£) , w 2 (£) , w 1 _(^), w 2+ (£) , v + (£), and v_ (£) are the counterparts 
of Y 1 (k), Y 2 (k), Y^fk), Y 2+ 0O • U + <k) , and p_(k), respectively, and are 
defined as: 


2 2 
w^(o = l-r 


w 2 (C) = (l-MC) 2 - c 2 


w <o = u+o 


W 2+ U) = 


2,nv_ (£) = - 2nM ( l+£ ) 1/?2 - i. J 


-1 arctan 


(1-mV) | Wl (v) 
| W 2 (V) 


V + (0 = 


|w 2 (C)+w 1 (^) ( 1 -M^) 2 }v_(?) 
(?-u(M)) (^-u* (M) ) 


The branch cuts of w^(£) lie from -1 to and from 1 to 00 and the branch 

cuts of w-(C) lie from r^-r to -oo and from rrrr to 00 in the £ -plane, which 
2 M— 1 M+l 

_ig 

correspond to the branch cuts of y (k) from -k to ooe and from k to 

1 k ° o 

coe^, and the branch cuts of Y- (k) from —7- to ooe ^ and from - — 77 to coe^ 

2 M— 1 1+M 

in the k-plane, where 


8 = arg k 

o 
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Since this rotation of the branch cuts does not affect the radiation 
conditions as discussed in Chapter 4, this is permissible. The contours 
of the integration C v and are shown in Figure 7-1. In the following 
figures we draw two contours on the same plane because after the inte- 
gration with respect to V, the pole V = £ is the contour of the next 
integration and therefore their relation can be easily shown. 


Let 


x = - (l-M^) r cos© 
o ' 1 o o 

y — — (l-M 2 ) * L//2 r sin0 
o ' 1 o o 

r and 0 are given as: 
o o r 



y o Vl -« 2 

0 = arctan 

o x 


r > 0 
o 


7T > 0 > 0 . 

O 


(7-9) 


(7-10) 



Figure 7-1 Branch Cuts and Integration Paths 
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Deform the contour into the curve Oo^ defined by 
cos (9 +iT)-M 

V = 2__ * (7-11) 

1-M 


where T is real and runs from -<» to oo. This curve is a branch of a 

±i9 0 

hyperbola with asymptotes V=e and has a vertex at 


V 


cos0 -M 
o 


1-M 


2 " 


(7-12) 


This vertex lies between the branch point V = -1 and ^ 


cos0 > M + M - 1 . 
o 


(7-13) 


Therefore, if Eq. (7-13) is satisifed, the contour can be deformed into 
the hyperbola without cutting the branch lines and we restrict the 
position of the line source in this region* The region expressed by 
Eq. (7-13) depends on the Mach number M, and as is shown in Figure 7-2, 
when M is small it covers almost all the moving flow and when M 
approaches 1 it is restricted to a small region. 



Figure 7-2 The Region Expressed by Eq. (7-13) 
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The contour of integration is deformed m a similar way. Let 
x = rcos0 ) 

(7-14) 

y = rsin0 ) 

with r >0 and *rr > 0 >0. Deform the contour into the curve 

£ = cos(0+iT^) (7-15) 


where is real and runs from -co to oo. This curve is a branch of a 

±i0 

hyperbola with asymptotes £ = e and has a vertex at 


£ = COS0 . 


(7-16) 


This vertex lies between the branch points £ = -1 and £ if 

1+M 


COS0 < 


1 

1+M # 


(7-17a) 


If 


COS0 < 


1 

1+M 


(7-17b) 


the contribution from the branch cut must be added. In both deformations, 
no contribution occurs from ‘the linking arcs at infinity. Then on the 
contours the exponential part of the integration (7-6) becomes 

-lk (vx +w^(v)y -Cx-w, (£)y ) = -ik f-r coshT+r Mcos0 -rcoshr ) . (7-18) 

o\o2 a o 1/ o\o o o 1/ 

The integration takes a different form according to the geometrical 
relation between the position of the line source and the observer as 
follows . 
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Figure 7-3 Region Expressed by Equation (7-21) 
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Region 1 


When 


COS0 < 


cos0 -M 
o 

1-M 2 


(7-19) 


Introducing the real geometric angle and length defined by: 


x = -r, COS0, 
o 11 

y o = -r «lnO 


(7-20) 


Eq. (7-19) can be written as: 


0 > 0 . 


(7-21) 


where 


COS0, 


COS0 2 = 


M 


2 2 2 
1-M sin 0 (1-M ) 


1-M 


(7-21a) 


The restriction of the angle (7-13) ensures that Eq. (7-21a) always has 
the root. The region expressed by Eq. (7-21) is shown in Figure 7-3. 

If Eq. (7-19) is satisifed 

cos0 > cos0 (7-22) 

o 

since 

COS0 -M 

COS0 > “ T“ (7-23) 

° 1-M 2 

is satisfied identically and hence 

9 < 9 . (7-24) 

o 
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AS 


cos9 — M 
o 


1-M 


1 

1+M 


is also satisfied identically 


cos0 < 


1+M 


(7-25) 


(7-26) 


and hence no contribution occurs from the branch cut. The relation of 
the two contours is as shown in Figure 7-4. 

The transmitted sound field <|> is given as: 

4^ (x,y ,w) = 4> d (x,y,a)) (7-27) 


where 




oo 

-f 


ik rcoshT. 

F (?)e ° 1 

2~2 2 dT l 

4tt a (C-u) (C-u*)w. (C)w, (C)v (C) 

1 + 2 + + 


(7-28) 


and 


. p o (« 


^ V_(v)w i 


V_(v)w i (Vjw^^fv) (l-MV)e 


ik r (coshT-Mcos0 ) 
o o ' o' 


2tt(V-£) 


Vi-m 2 


-dT . (7-29) 


If and intersect as shown in Figure 7-4 (b) , which may occur when 
0 ~ 0 2 , the contribution from the pole V = £ must be added to give: 

(x,y ,(a)) - ^(Xry^w) + <() (x,y ,u>) (7-30) 
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(b) 


Figure 7-4 


Relation Between C 


a 


and C V1 


0>e 2 
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where 




(l-MOw^le 


, 2 2 
4ir a 


./ 


(l-MC)w (€)e 


ik (rcoshx - x -w„(C)y \ 

O' 1 O 2 O' 

|w 2 <C)+w 1 (C) (l-MC) 2 j- 1 

ik (rcoshT -Cx -w _(C)y \ 

O' X O 2 Of 


A 2 2 
4ir a 


|w 2 (?)4w 1 (C) (1-MC) 2 | 


dT 


1 • 


(7-30a) 


X , X 2 , X 3 , and X^ are the roots of 


cos(9 +iT)-M 


cosCQ+it^) 


and 


1-M 


\ > X 2 > X 3 > h * 


(7-30b) 


(7-30c) 


Let 


cosCS^+iT^) = u(M) . 


IT 


(7-31) 


As shown in Appendix A, 0 3 =— . If 0<0 3 , the contribution from the poles 
£ = u(M) and C=u*(M) must be added to give 


= <J , d +H{ V 0)<{, pi 


(7-32) 


where 


^Pl 


F (?)e 
o 


ik Q (Sx+w^Oy) 


2ira (C-u*) W;l _(?)w 2+ (Cjv + (^) 


C=u 


P o (C)e 


ik Q x+w^ K)y) 


2Tra (C-u)w 1 _(C)w 2+ (Ov + (?) 


(7-32a) 


S=u* . 


The relation between the poles and the contours is shown in Figure 7-5. 
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Figure 7-5 Relation Between the Poles 
and Contours 


Region 2 
When 

0 < 0 . (7-33) 

o 

The region expressed by Eq. (7-33) is shown in Figure 7-6. If Eq. 
(7-33) is satisfied, then 

cos0 -M 

cos0 > 2_ . (7-34) 

1-VT 

The relation between the two contours is as shown in Figure 7-7. When 
we deform the contour to the two contours intersect and the pole 

£=V brings another sound field (f> t which is the same as the solution for 
the doubly infinite vortex sheet problem. In the calculation of we 
are free to deform the contour if there is no contribution from the 
linking arcs at infinity as long as lies in the left of C^. We 
calculate <f> along the stationary path £(a) which has its saddle point 
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Vortex Sheet 


Source 




Figure 7-6 Region Expressed by Equation (7-33) 



Figure 7-7 


Relation Between C and C 9 < 0 
VI £1 o 


i ±9>i 

between -1 and t— - on the real axis and has its asymptotes £ = e q 


1+M 


where : 


9. = arctan 

4 


-(l-M 2 ) 1/2 y +y 
o 

x-x 

o 


(7-35) 


The character of £ (CT) is discussed in detail in [l]. Then 


♦t- 


00 

/ 


i (1-M£)e 


ik g (a) 


4ir 2 a 2 |w 2 (?)+w 1 (?)(1-mC) 2 1 


d£ 

da 


da 


where 


g(a) = -£x -W (£)y + £x+w (£)y Img(O) =0 . 

O Z O 1 


(7-36) 


(7-36a) 
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If C (d) captures the poles u(M) and u* (M) , another sound field <j>p 2 roust 
be added. 



(l-MC)e 


2iTa 


(1-M5)e 

„ 2 
2TTa 


ik (?x+w (£)y-?x -w (?)y ) 

O \ 1 O 2 O / 

^■|w 2 (?)+w 1 (?) (1-MQ 2 J 


?=u 


C=u* . 


(7-37) 


does not appear when 9 q < 9^ where 


cos(9 5 +H 5 )-M 


1-M 


2 


= u(M) . 


(7-38) 


The relation between the contours and the poles is as shown xn Figure 
7-8. If 0 < 9- , where 

b 

cos6 ff = (1+M) ~ 1 , (7-39) 

b 

the contribution from the branch cut d>. , must be added. 

bl 




the Poles and the Contours 
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(7-40) 


ik rcoshT, 

F «)e ° 1 

— 2 dT l- 

47T a « -u) (C -u* ) w 1+ (C ) w 2+ (C ) v + « ) 

[f(C)l = f (C+ie) -f(£-ie) is the discontinuity of f(S) across the 
branch cut. 

As a result, the transmitted sound field <|> in Region 2 is given as 

^l = ^d +H (l" 9 )^Pl +H(0 6" e) ^bl + \ +2(M ' x ' y)H(6 o" 0 5 ) ^P2 ' (7_41) 

where Q(M,x,y) is 1 when C(a) captures the poles and otherwise zero. 

<J>h* "$p 2 are the solutions for the doubly infinite vortex sheet problem 
(without rigid plate) . 

Region 3 
When 

9 > 0 > 0 . (7-42) 

2 o 

Region 3 is between Region 1 and 2, and is shown m Figure 7-9. 

The relation between the two contours is as shown m Figure 7-10. The 
contribution from the pole <|> partially appears. 

dT 1 . (7-43) 


K tl 


/ 


(1-M^)w 1 (?)e °' 


lk (rcoshT, -?x -w 


o 2 


4ir 2 a 2 |w 2 (?)+w 1 (?) (1-MC)‘ 


*bl 


-f 


x(8-9 6 ) 
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0.8 



Figure 7-9 Region Expressed by Equation (7-42) 



Figure 7-10 Relation Between C and C 

0 > 6 > 0 V1 w 

1 o 


X and X are the roots of Eq. (7-30b) , and 
o o 


X 6 >X 5 


x e - - x 5 • 


(7-43a) 


If 0 < 0^ , the contribution from the branch cut ~ must be added. 
6 2 


b2 


-/ 


l( e-e 6 ) 


i 




ik rcoshT -£x -w_(£)y 

~ 1 1 O 2 


°\ 


„ 2 2 
4ir a 


|w 2 (C)+w x (^) (1-MC) 2 | 


dT 1 . (7-44) 


If 0 < - and 0 q > 0^_ f the contribution from the poles <|> appears. As a 
result, the transmitted sound field in Region 3 is given as: 

*1 = + H (t" e )*Pl +H(6 6- 9, *bX + *U +H(6 6- e, *b2 + H (|- j « (8 o- 0 5>*P2 


(7-45) 
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CHAPTER 8 


ASYMPTOTIC EVALUATION OF FAR FIELD 

The fields <j)., 6 , ({>_ , <f>, . , and (f), 0 can be evaluated asymptotically 
d t tl Dl DZ 

under the assumption that k r » 1 and r » r . First consider the case 

o o 

when 0>0 . It is permissible to change the order of integration in 
o 

integrating (7-28) and (7-29) because and T are real. Now assuming 
k Q is real and positive since T^ = 0 is the saddle point, a standard 
stationary phase method is used to evaluate the integration to give: 

ik r+yi 
F (Oe ° 4 

<t> = =75-5 ° — (8-1) 

a 2-/1 it j/ a (S-u) (C-u*)w_ (£)w-(Qv(i;)A r 

1+ 2+ + O r n 

4=COS0 . 

If 0 < 0 _ , C-. must be deformed to enclose the branch cut as shown m 
6 £1 

Figure 8-1, and special care must be taken since w 2+ ^' w 2 ^' 

take the different values depending on from which side C approaches to 

the branch cut. 
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V*bi 


lk rcoshT 
F (C)e ° 1 

—22 2 dT l 

47T z a z (C-u) (S-u*)w_ (Qw 0 (Qv (C) 

1 + 2 + + 


ik rcoshT\ 
o 1 


F (C)e 
o 


' 4ir a (?-u) (C-u*)w 1+ (C)w 2+ (C)v + (C) 


aT 1 + o(|k o r |- 3/2 ) 


F (?)e 
o 


ik r + — i 
o 4 


4/2 ir 3/2 a 2 (?-u) (?-u*)w (?) vic~r 

1+ o 


£=COS0 


w 9 (C)v (C) I W (C)v (C) i 

C=-O+cos0 C--K)+cos0 


ta- 2 ) 


where for example, f(£) | means the value of f (£) when £ approaches 

C=+0+a o 

a Q through positive value on the upper side of the branch cut. As 0 = 0 
is the saddle point of integration (7-36), we have: 


ik g(a) +sgng" (a) “i 

j m i(l-MOe ° 4 

2/2 ir 3/2 a 2 |w 2 (C)+w 1 (?) (1-M?) 2 j. /iTJ? 


(8-3) 


Since we assume r » r , r =0 is the saddle point of integration (7-43) 

o 1 

and the mam contribution to the integration comes from this point. 
Therefore, asymptotically it is permissible to deform to £(a). We 


*ti + *b 2 
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It is concluded that Region 2 and Region 3 asymptotically have the same 
character . 

As we assume r » r , both C„. and £(a) have the same asymptotes 
o £1 

and for small Mach number the poles are not included in the deformation 

7T 

from to £(0) if 9 is not so close to Therefore, we have: 

\ + 5(M.*.y>H(0 o -9 5 )? p2 = + a(j-e)H(e o -0 5 )* p2 

- + +b2 +H (T- 0 ) B(6 o- e 5^P2- <8 - 5) 
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CHAPTER 9 


THE IMPULSIVE PROBLEM 


In this chapter inverse transform with respect to go is carried out 
to obtain the solution for the impulsive problem, ^(x^t) is formally 
given as: 


<f> x (x,y,t) 



-T- , v -I0)t_ 

(}> 1 (x,y,oo)e doo 


(9-1) 


where integration is carried out along the real axis m the co-plane. 
In Region 1 we have 


= ^ + h (t _ e ) ^ 


pi 


(9-2) 


V 


uu w 

■// 


v (V) w (V) w (V) (l-MV)S<t-(r coshT-r Mcos0 +rcoshT n ) /a> 
1- 2- ( \ o o o 1/ ) 


4tr 2 a 2 Vl-M 2 (V-g (C-u) <C-u*) Wl (Ow^lOv, (C) 


dTdT, 


1+ 


2 + 


(9-3) 


where the integration with respect to co has been performed. This can be 

integrated with respect to either or T, but the remaining integration 

cannot be performed explicitly. It is however clear that (J>_ vanishes 

d 

for times at < r - r Mcos0 + r. The arrival time of the wave front 
o o o 

becomes earlier or later depending on the position of the line source 

due to the term r Mcos0 and hence the existence of the flow in the 
o o 

lower half-plane. 
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dxdto 


^P1 




v (v)w, (v)w 2 _(v)(l-Mv)e 




coshT-Mr cos0 +ux+w. 
o o 


1 (u)y|/a| 


4ir 2 a 2 (v-u) Vl-M 2 (u-u*)w (u)w (u)v + (u) 


V (V) w, (VJw^ 




iu) f-( r o 


coshT-Mr cos0 +u*x+w 
o o 


1 (u*)yj /aj. 


-dXdU) 


4ir 2 a 2 (V-u*) Vl -M ^ (u*-u)w, (u*)w_,(u*)v (u*) 

1 — 2 + + 


(9-4) 


It can be expressed in terms of the ultradistribution S(z), where z is 
complex, as: 


d> = 
y Pl 


CO 

/ 

_oo 

00 

/ 


v_(V)w (V)w 2 _(V) (1 


-MV) <5 {t-( r Q 


coshT-Mr cos0 +ux+w. 
o o 


1 (u) y) /aj- 


2 / 2 

27Ta (V-u) Vl-M (u-u*)w, (u)w_, (u)v t (u) 


dT 


1- ' 2+ +' 


v (V) w_ (V) w (V) (l-Mv)6<t-(r coshT-Mr cos0 +u*x+w 
- 1- 2- ( \ o O O 


1 (u * ) y)/a} jT 


27ia^ (V-u*) Vl-M^ (u*-u)w (u*)w (u*) v (u*) 

1— 2+ T 


(9-5) 


The ultradistrxbution 5 (z) has the property that 


6 (z) = E (i Im z) n 6 (n) (Pez)/n: (9-6) 

— OO 

and since generalized function 6 ^ (Rez) vanishes for Rez ^ 0 , indi- 

cates the instability waves and vanishes for 

at < r - Mr cos0 + Reux + Rew_ (u)y . (9-7) 

o o o 1 
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It is easily checked that Reu f Reu*, Rew^fu), and Rew^fu*) are all 
positive and these instability waves appear only in the downstream. 
Therefore (f)^ satisfies causality. On the other hand, if we deform 
the contour m the k-plane so that the two poles are always below the 
contour as discussed m Chapter 4, <J) pi appears in the upstream and for 
the certain value of negative x, <f) pl can be heard t<0 and causality is 
violated. 

In Region 2 we have 

<t>l = < f , d +H Gf- 0 ) < J ) p 1 + H( V 6)<f bl + ♦t + Q( M ' x «y)H( 0 o - 0 5 )4, P2 (9 " 8) 


00 1 ( 0 - 0 .) 

W/ 

-00 0 


v (Vjw^ (V)w 2 _(V)(l-MV) filt-^r^coshT-Mr^cosS^+rcoshT^j/al 
4Tr 2 a 2 (v-^) Vl-M 2 (C-u) (?-u*)w 1 .(C)w_. (5)v x (?) 

• 1+ Z+ + 


dT x dT. 


(9-9) 

6 . can be integrated once only, but it is clear that it vanishes for 

bl 

times 


at < r - Mr cos0 + rcos(0 -0) . 
000 6 


(9-10) 


Since it is restricted to the region 0 < 0 , it is also causal. As 

6 

indicated m Chapter 7, (j)^ and <J) are the solution for the doubly 
infinite vortex sheet problem (without the rigid plate) which is calcu- 
lated exactly and discussed in detail in [6]. The only difference is 

that the instability waves <j> _ are blocked and do not appear when 9 <0„. 

P2 o 5 

In Region 3 the integration cannot be performed explicitly for <j> and 

<j) However, as shown m Chapter 8, its nature is asymptotically the 
b2 

same as that of Region 2. 
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CHAPTER 10 


CONCLUSIONS 


The transmitted sound field satisfies the full Kutta condition 
and causality condition. 

The sound waves (f>^, (f>p^r and cf)^ can be heard first at the edge of 
the plate at the time 

t = (r -Mr cos0 )/a (10-1) 

o o o 

which is the time that the emitted sound wave takes to travel from the 
source point to the plate edge in the moving fluid. This arrival time 
becomes earlier or later depending on whether the source is in the 
upstream ^0^ < —j or m the downstream of the plate edge 

The field 0^ represents the field scattered from the edge of the 
plate and can be heard everywhere in y>0. The wave front clearly 
occurs on the circle 

r = at - r +Mr cos0 (10-2) 

o o o 

and the time t = (r+r -Mr cos0 )/a is the time for a signal to travel, 

o o o 

via the edge, from the source to the observer. 

d>. . is the bow wave associated with d>, and can be heard in the 
bl T d 

region 0 < 0 . Its wave front is expressed as : 
o 

at = r - r cos0 + rcos (0 -0) . (10-3) 

o o o 6 
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This is the straight line from x = (1+M) (at - r + Mr cos0 ), y = 0, to 

o o o 

the point of contact with the cylindrical diffracted wave front. <f> p ^ 

can be regarded as an instability wave triggered by the diffracted wave 

7T 

at the edge and is restricted to the region 0 < — . It exists m the 

triangular region bounded by the line x = y and Reux + Rew^ (u)y = 

at - r + Mr cos0 . 
o o o 

The transmitted sound fields can be classified depending on how 
<{>t wave is transmitted. <j> t and <|> are the transmitted waves through 
the vortex sheet and are identical m the form to those found for the 
infinite vortex sheet. 

In Region 1, <|> is completely shaded by the rigid plate. In 
Region 2, <f> is fully transmitted. Region 3 is the transition region 

between the two and (j>^ is partially transmitted depending on its wave 

number. 

If the observer is close to Region 1, <J>^ is mostly shaded and if 
the observer is close to Region 2, <J> t is mostly transmitted. If the 

observer is far from the plate edge, the nature of the sound field is 

the same as that of Region 2. 

<{)p 2 is another instability wave which occurs when £ (a) captures 
the poles. If the source is close to the rigid plate (0 < 0^) , this 

wave is blocked by the plate and does not appear. 

If M tends to zero, the boundaries of Regions 1, 2, and 3 approach 

the diffracted angle of the Sommerfeld problem 
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APPENDIX A 


POLES OF K 


The poles of K are given as the roots of the equation 
Y 2 (k) + y 1 (k)a 2 = 0 . (A-l) 

Introducing the parameter 

S = ka/£ , (A-2) 

Eq. (A-l) reduces to the form: 

+w 2 (?)| +w i (?)w 2 (?)J = 0 (A-3) 

where 

wJ(C) = 1-C 2 

2 2 (A " 4) 

w 2 (?) = (1 -mo -r 

and to satisfy the sign convention of and y^, the imaginary part of 
their counterparts and w 2 must have the same sign as k Q along the 
real axis. Accordingly , the branch cuts in the 3-plane are as shown 
in Figure A-l. 

i) w (C) + w 2 (£) = 0 (A-5) 

V^(C) = w 2 (0 

2 2 2 

l - c = (l-MC) - Z 

(1 - 1 + M3) (1 +1 - M3) = 0 
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3-plane 



© 



o 

□ 


indie ate s the 
arg. of A -3 

indicate s th e 
arg of /i+3 


Figure A-l(a) 


Branch Cuts of w^(3) 



C = 0 



It is apparent that 
2 

L - — is situated on 
M 


£ = 0 is not the solution of Eq. (A-5). If M<2, 
the right of £ = 1 and cannot be the solution of 


Eq. (A-5) . 


(A-6) 


11) C 2 + w 1 (?)w 2 (C) = 0 

4 2 2 

C = wjiciwjic) 

1 - 2M£ + (M-2) C 2 + 2M? 3 - M 2 £ 4 = 0 




= 0 


Z = 


c " 


-M i V M 2 +4+4 Vm 2 +1 
2+2 VM 2 +1 


M ± Vm 2 +4-4 Vm 2 


+1 


: Vr 2 


2 V 1+M -2 


S ■ 


-M+ V M 2 +4+4 V 1+M 2 

:V?~ 


is between £ = 0 and £ = 


1+M’ 


c- 


-M-V 


M 2 +4+4Vm 2 +1 


2+2 VM"+1 2+2>/l+M 2 

is between C=0 and £=-l. These two are not the solution of Eq. (A-6). 


C = 


±Vm 2 +4-4 Vm 2 +1 
2 Vl+M 2 -2 

V j - 2 


are the complex numbers which are situated above 


and below £ = 
of Eq, (A-6) . 


1+M +1 r - 

, respectively, when M < 2/2 and are the solution 


2M 


These two solutions can be written as: 


C - 

where T is 
coshT 


COS 7T + irj 

the positive root of the equation 
_ 1 +Vi+m 2 

1/2 M 


(A-7) 


(A-8) 
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APPENDIX B 


CALCULATION OF y (k) AND y_ (k) 


As y (k) does not have any zeros , y (k) is defined by 


In y (k) 


-1 f £n y (X) ^ 

= 2ttT J r-k ' dA 


(B-l) 


where C^ passes below X = k. Introducing the new parameters 


? ■ 

v = XA 0 , 


(B-2) 


Eq. (B-l) can be written as 


£n u (k) 


- f 

2th J 


Jin p(Vk ) 

o__ 

V-? 


dV . 


(B-3) 


The branch cuts of w (V) lie from -1 to -°° and from 1 to 00 and the 

branch cuts of w- (V) lie from -~r to - 00 and from r^TT to 00 in the V-plane 
2 M-l M+l * 

as shown in Figure B-l where w (V) and w (V) are the counterparts of 

-L ^ 

y^(k) and y 2 (k) and are defined by: 


W 2 (V) = 1 - v 2 

W 2 (v) = (1-MV) 2 -V 2 . 


(B— 4) 


The contour C^ is distorted so that it is wrapped around the cut from 
-1 to - 00 . This process gives an additional contribution from the 
semicircle of radius R. 
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Figure B-l Branch Cuts and Integration Paths 


As | V J ->°° in the upper half-plane, w =-iV + 0(l), w =0(|v|) and so 

I * 


y(k) 


-iM 2 y 

k 

O 


+ 0 ( 1 ) • 


(] 


In the lower half-plane w =iv+0(l), = 0 ( | V | ) and so 


y(h) = 


. 2 
iM V 

k 

o 


+ 0 ( 1 ) . 


(1 


It follows that the contribution from the semicircle is: 



-*/«■( 


. 2 i0 N 
•lM Re 


1+0 (1/R) 




1+0 (l/R) d0 


I +0 (^) 




The discontinuity in Any (k) across the cut is : 


w 2 (v+ie) +w 1 (v+ie) a 

. - 

w 2 (v-ie) +w^ (v-ie) a 


An(-l) 


(l-M V j |w 1 |i+|w 2 | 

An r— — = 2 i arctan 

-(l-M V j )w 1 |i+|w 2 


/, 2 2 \ 
(l-M V ) 


•< V < -1 . 


Zn\x (k) can be given as: 


-1 arctan 


(i-mV] 


„ „ , 1. M R 1 C dV 1 C 

An U _(k) =- I An— -- y — ? ^ 


— 7 to T 5 -? to< - l -0 + 7 to( -*> -7 / 


-1 arctan 


(l-M 2 V 2 ) [wj L 
| W 2 1 “ 2 
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= - In 


M(l+C) 1/2 

k 1 / 2 


1 

7T 


-1 



arctan 


(l-M‘ 

> 9 

\T 

ll w ll * 


W 2 

2 


v-C 


dV . 


(B-9) 


The asymptotic behavior of y_(k) as |k| -*■“ is easily found as 

y_(k) = 0 (|k|- 1/2 ) . (B-10) 

The identity y + (k) = p(k)y (k) can then be used to derive the corre- 
sponding results for y + (k) as follows: 


V + 0O 


(y 2 +Y 1 a 2 )u_(k) 

(k - u(M)k Q j (k-u* (M)k_J 


(B-ll) 


The asymptotic behavior of y + (k) is 
y + (k) = 0 ( | k | 1/2 ) . 


(B-12) 
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APPENDIX C 


SOLUTION FOR THE SOMMERFELD HALF-PLANE DIFFRACTION PROBLEM 


Consider the case m which there is no flow in y<0. The 
governing equations are : 

9 2 <J> 

— ~ - a 2 V 2 (j) = 6 (x - x ) 6 (y - y ) 6 (t) y < 0 

9t 2 2 ° ° 

- aV<j> =0 y > 0 . 

9t • 

Carrying out the Fourier transform defined by Eq. (2-2) , Eq. (C-l) 
reduces to the form: 



3 <j>„ 


3y 


+ Y (j) 

2 y r 2 


6 (y-y ) -ikx 
o o 

2 2 e 
4^T a 


3y' 


- + Y 2 ? 

2 'll 


= 0 


where 


2 , 2,2 

Y l*V k , 


k = 03 /a 
o 


j (C-2) 


Solutions to the set of ordinary differential Eq. (C-2) can be easily 
found as: 
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_ iY,y -iY,y 

4>1 = A^k^Je +B 1 (k,w)e 


= -iY,y -iY,y 

4> 2 = A 2 (k,u)e + B 2 (k,w)e 



The sign convention of ** s chosen such that the sign of its imaginary 
part is always positive and therefore and B 2 must be set to zero. 

The boundary conditions at y=0 are: 


3 ^ 3*2 

cStT = 3t“ 

3 ^ 34>2 

3t“ = dt~ 


34> x 

3y 



x>0 (continuity of pressure) 
x > 0 (continuity of displacement) 
x<0 (rigid plate condition) . 


(C-4a) 

(C-4b) 

(C-4c) 


To express these conditions in the transformed region, carrying out 
the half-range Fourier transform defined by Eq. (2-7) we obtain: 




3^. 

9y 



3<f> 2 _ 

3y 9y 0 
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Using the relation 


w ■ l Vi 


^2 = e l(kX o +Y i y 0 t 

ST ■ “’A • 4 *V 


and eliminating all the unknown plus functions except for 


obtain : 




- 1 (kx +Y y) 

O 1 

5 

m 

4tt a 


(C-7) 


The branch cuts of y^ are the same as those defined in Chapter 3. The 
common region of analyticity is 


link > link > - link , 
o o 


and Ima) is kept to be positive. The function y^ is split into two 


parts such that y^ = y^ /y^ + where 


y i + “ (v k )' 


<!- - ( k o +k ) 1 


Plus functions are analytic in the region link < Imk^ and minus functions 
are analytic m the region Imk>-Xmk o . The Wiener-Hopf Eq. (C-7) can 
be rewritten as: 


ly 1 = 1 V k o +k (♦l--*2-) ' 72^ 


i (kx +y.y ) 
o ' 1 o 


4tt a yk^-k 
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The last term of Eq. (C- 10) can be decomposed into the sum of plus 
and minus functions as: 


-i (kx +Y y ) 
o 1 o 

F (k) = - — F (k) +F i_ (k) ' (C-lla) 

4TT a 


where F^ +/ F^ are given as: 


f 

2TTi J 

C+ 


F i+ (k > -2 ir / tI- 


? i- (k > ■ hi f 


F(X) 

X-k 


dX 


c- 


(C-llb) 


The contours C and C lie in the common region (C-8) and C, passes 
+ — + 

above X=k and C passes below X = k. The exponential decays of F ( X) 

on the contour/ since Imy., >0, y < 0 r ensures that these integral 

1 o 

exist and that F^ + (k) = 0 ( | k | (k) = 0^|k| 1 j as k*^°° in respective 

regions. 


The Wiener-Hopf equation now takes the form as: 


2 



9y 


+ F 


i+ 




F x _ = C(k) . 


(C-12) 


The function C(k) is the entire function of k and must be a regular 
function of k in the whole k-plane. At the plate edge, we assume: 
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as x +0 on y = +0 


(C-13a) 



<f>i C ± (w) 
^2 ^ C 2 (a)) 


as x -0 on y = +0 
as x -► -0 on y = -0 • 


(C-13b) 

(C-13c) 


Carrying out the half-range Fourier transform defined by Eq. (2-7) , 
the asymptotic behavior of the above functions can be calculated with 
the aid of the Abelian theorem as: 


*1 + n( 

1, | -l/2\ 




3y ( 

|k| ) 

as 

k 

->■ 00 

= / 

i l~l\ 





|k| M 

as 

k 

-> 00 

1 

II 

O 

I*!" 1 ) 

as 

k 

-»■ 00 

Hence from Eq. 

(C-12) 




C(k) = o | 

W 1 ) 

as 

k 

oo 

C(k) = o| 

i*r 1/2 ) 

as 

k 

-*■ 00 


m 

Imk 

< Imk 

o 

(C-14a) 

in 

Inik 

> -Imk 

o 

(C-14b) 

m 

Imk 

> -Imk 

o 

(C-14c) 


in 

Imk < Imk 

o 

(C-15a) 

in 

Imk > -Imk 

o 

(C-15b) 


C(k) is regular m the whole k-plane and tends to zero as k tends to 
infinity in any direction. Hence from Liouville's theorem C(k) must be 
identically zero. Therefore, we obtain: 


3<t>l + _ -Vk 0 -k VO) 
3y 2 


(C-16) 
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A 1 (k r 0)) = (J) 1 j 

y=0 



3y 


iF i^ (k> 

YV 1 


(C-17) 
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APPENDIX D 


THE SOLUTION FOR THE INFINITE VORTEX SHEET PROBLEM 


The governing differential equations are the same as Eq. (2-1) . 

The boundary conditions (2-4) , (2-5a) , and (2-5b) must be satisfied 
throughout the vortex sheet. Instead of the half-range Fourier trans- 
form, using the full-range transform, Eqs. (2-4) , (2-5a) , and (2-5b) 
can be written as : 

a? 2 = (D-D 

3 ?, 

— = -laxxh (D-2) 

dy 

% 

- 5 — = -iioh . (D-3) 


From Eqs. (D-2) and (D-3) we obtain: 


a 


3<J) 2 

8 y 8 y 


From Eq. (2-12) 


= A, 


y =0 


-xkx 

e sinY_y 

= A 2 + 2~2 ~ 

y=0 4iT a Y 0 


(D-4) 
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y=0 


^2 

3y 

y=0 


- 1Y 2 A 2 


-ikx 

e °cosy y 

2 O 


„ 2 2 
4tt a 


(D-5) 


Substituting Eq . (D-5) into Eqs. (D-l) and (D-4) gives 


-ikx 

a e °sinY 9 y 
4TT a Y„ 


-ikx 

o 

e cosY y 

lCtY i a i - - lY 2 R 2 -22 

4TT a 


(D-6) 


Finally, from Eq. (D-6) , and A^ are determined as: 


\ = 


ia e 


-ikx -iY,y„ 

O 2 o 


2 2 2 
4tt a (Y 2 -kx y x ) 


(D-7) 


-ikx 


A 2 = 


i e 


0 2 2 
87T a y 


A Vo Y 2‘ aY l " iY 2 y c 


Y 2 + a Yj_ 
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